Double-layered monopolar order in Tb2Ti 2 07 spin liquid 
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Ho2Ti207 and Dy2Ti207 spin ices exhibit elementary excitations akin to magnetic monopoles. 
Here we focus on Tb2Ti2 07 spin liquid, where correlated magnetic moments keep fluctuating down 
to very low temperatures. Using a monopole picture, we have re-analyzed the field-induced magnetic 
structure previously determined by neutron diffraction in Tb2Ti207- We show that under a high field 
applied along a [110] direction, Tb2Ti207 orders as a three dimensional arrangement of monopole 
and antimonopole double layers. In contrast, Ho2Ti2 07 spin ice in the same conditions behaves 
as a monopole-free state. By symmetry analysis we derived the distortions compatible with the 
observed magnetic structure of Tb2Ti207 which can be related to the appearance of the double- 
layered monopolar order. 

PACS numbers: 71.27.H-a, 61.05.fm, 75.25.-j 



I. INTRODUCTION 

Spin ices are geometrically frustrated magnets, which 
support exotic ground states and excitations^. Their 
zero-field magnetic ground state has an extensive de- 
generacy issued from topological constraints, like those 
which govern the two energetically equivalent positions of 
protons of water molecules in real ice. Spin-ice behavior 
has been found in the pyrochlore compounds Ho2Ti2 07 
and Dy2Ti207, where the rare-earth (R) ions reside on 
a lattice of corner sharing tetrahedra, and the magnetic 
moments behave as Ising spins, being constrained to lie 
along their local (111) anisotropy axes 2 -* 3 -. The R mag- 
netic moments interact via dipolar and superexchange 
interactions, yielding an effective nearest-neighbor (nn) 
interaction of ferromagnetic (FM) nature which governs 
most of the physics. In spin ices, the local structure which 
minimizes the energy of a single tetrahedron obeys the 
so-called ice rules, having two magnetic moments point- 
ing inwards and two outwards of the tetrahedron (Fig. [TJ 
left panel) . This situation gives zero net magnetic charge 
on the center of each tetrahedron. Starting from this lo- 
cal "2-in, 2-out" structure, emergent excitations can be 
created by flipping one magnetic moment (Fig. [TJ right 
panel). These point defects with "l-in, 3-out"/"3-in, 1- 
out" configuration of one tetrahedron can be considered 
as magnetic charges 4 or monopoles 5 . A spin flip induces 
two oppositely charged monopoles connected by a Dirac 
string 5 - when taken apart. Because of the extensive de- 
generacy of the spin-ice ground state, it takes only a finite 
energy to separate these monopoles to infinity 5 . The en- 
ergy cost necessary to violate the ice rules and create a 
monopole pair can be provided either by thermal fluctu- 
ations or by applying a magnetic field 6 . 

Tb2Ti2 07 spin liquid has a different zero-field ground 
state, akin to a cooperative paramagnet where strongly 
correlated magnetic moments fluctuate at short time 
scales down to very low temperatures (50mK) 7 . The 



exact reason why it behaves as a spin liquid is still a 
matter of debate. Virtual crystal field excitations have 
been previously invokec&i£. Alternatively, we recently 
proposed to take into account the influence of a distor- 
tion, which induces a two singlet ground stated. In any 
case, two key differences between Tb2Ti207 and Ho/Dy 
spin ices must be outlined: (i) The first excited crystal 
field doublet lies at energy above the ground state dou- 
blet which is one order of magnitude smaller than in the 
spin ices, giving the magnetic moments a certain degree 
of freedom with respect to the local (111) axes; (ii) The 
balance between exchange and dipolar energies turns the 
effective nn interaction to be antiferromagnetic (AFM). 

In the present paper we study the long-range ordered 
magnetic structure which is induced in Tb2Ti2 07 un- 
der a field H applied along [110]. We have re- analyzed 
the magnetic structure previously determined by neutron 
diffractio n 12 ! 13 using the monopole description. We show 
that above a critical field the latter magnetic structure 
which violates the ice rules can be viewed as a three- 
dimensional (3D) periodic arrangement of monopole ("1- 
in, 3-out") and antimonopole ("3-in, 1-out") double lay- 
ers. This structure which is unique for a given field re- 




FIG. 1. (Color online) Magnetic monopoles in spin ice. Left 
panel: Spins of two adjacent tetrahedra obeying the "two-in, 
two-out" ice rule. Right panel: Monopole-antimonopole pair 
induced by flipping of the spin connecting the tetrahedra. 
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lieves the spin-liquid degeneracy. In contrast, the spin ice 
Hc»2Ti207 in the same conditions selects a "2-in, 2-out" 
non-degenerate ground state which complies with the ice 
rules and can be considered as a "vacuum" state, namely 
a state without monopoles. Such description has interest- 
ing consequences when considering the spin excitations in 
the two compounds. A single spin flip creates either a va- 
cancy pair (in the Tb case) or a monopole pair (in the Ho 
case) . Whereas dipolar interactions seem to be sufficient 
to describe the monopole- free state of Ek^T^Oy spin 
ice 14 , the double-layered monopolar state of Tb2Ti2 07 
spin liquid cannot be explained by this mechanism. We 
propose that it could be stabilized by a distortion, either 
of Jahn- Teller or spin-Peierls type. 



II. EXPERIMENT 

Single crystals of Tb2Ti207 and Ho2Ti207 were grown 
by the floating-zone technique (see, e.g., Ref . [TBI) . Neu- 
tron diffraction studies were performed on the diffrac- 
tometer 5C1 (A = 0.845 A) at the Orphee reactor of the 
Laboratoire Leon Brillouin, Saclay. We used unpolar- 
ized neutrons, collecting typically from 200 to 400 re- 
flections up to sin#/A « 0.55 A -1 for each data set at 
T = 1.6K under the field H < 7T. The nuclear struc- 
ture parameters were deduced from the low-temperature 
measurements in zero field. The integrated intensities of 
these reflections were used to refine the components of 
the Tb 3+ /Ho 3+ magnetic moments using the symmetry 
constraints described below. The program FullProf 16 
was used to refine the nuclear and magnetic structure. 



III. MAGNETIC STRUCTURE ANALYSIS 

In rare-earth pyrochlores, the magnetic structure in- 
duced by a field H || [110] splits into orthogonal R-a and 
R-/3 chains, oriented along [110] and [110] respectively. R 
ions on the a chains have their local anisotropy axis at 
±35.26° from the field, whereas in the j3 chains the local 
anisotropy axes are perpendicular to H. In Tb2Ti207, 
when the field increases along [110], the ordering of the 
Tb moments occurs in two step G 13 i 17 i 18 . In the low-field 
region, magnetic order with k = propagation vector is 
induced. The magnetic moments in the FM-like a chains 
increase linearly with H and tend to saturate 6 /is) in 
a 1 T field. The FM-like order in the a chains means that 
the chains have non-zero homogeneous magnetic compo- 
nent parallel to the field. The moments induced on the /3 
chains remain quite small since the field is applied along 
the local hard axes of the Tb moments. Above a critical 
field (H c ~ 2T at T = 50 mK), the AFM order both 
inside and between the /3 chains occur. This is shown by 
the appearance of an additional k = (0,0,1) magnetic 
structure which yields magnetic intensity on the Bragg 
peaks of the simple cubic lattice (Fig. [2j middle panel), 
whereas below 2 T the magnetic intensity is restricted to 



the peaks of the face-centered cubic lattice (Fig. [2j left 
panel), in agreement with the Fd3m space group of the 
pyrochlore structure. The k = (0,0,1) structure cru- 
cially depends on the temperature T, the field H and its 
precise orientation A(/> with respect to the [110] axi s 13 i 18 . 

To describe these structures, we previously performed 
a symmetry analysis based on the theory of representa- 
tions and found that both k = and k = (0, 0, 1) mag- 
netic structures can be described by single irreducible 
representation s 13 ! 19 . This reduces the number of the re- 
fined parameters from 12 in an unconstrained refinement 
(3 for each of the 4 Tb/Ho moments in a tetrahedron), to 
6 and 4 in case of the k = and k = (0, 0, 1) structures, 
respectively. A detailed description of the analysis can 
be found in Ref. [l9|. 

An example of a typical fit for both k = and 
k = (0,0,1) structures of Tb 2 Ti 2 7 at 1.6K and 7T 
is shown in Fig. [3l The resulting magnetic structure of 
Tb2Ti2C>7 at 1.6 K in 7T calculated as a superposition 
of the k = and (0,0,1) structures is shown in Fig. |4j 
On each tetrahedron, both /? spins point either "in" or 
"out" (AFM behavior) and the a chain is made of al- 
ternating moments (FM-like behavior), so that the local 
spin configuration ("1-in, 3-out"/"3-in, 1-out") violates 
the ice rules. 

The high- field magnetic structure of Tb2Ti207 de- 
duced from this symmetry analysis agrees with our ear- 
lier results 12 and is consistent with the schematic picture 
proposed in Ref. [l8| (see Fig. 1) for Ising spins. We no- 
tice however that the model proposed in Ref. |l8| gives a 
magnetization value about three times smaller than that 
reported in the literature 2 ^. This discrepancy is likely 
to be attributed to strong extinction corrections which 
probably affected the intensities of 3 Bragg reflections 
measured with cold incident neutrons (5 A)A£. Our data 
analysis involves 200 to 400 Bragg reflections measured 
with hot neutrons (0.845 A) and weakly influenced by 
extinction. This lifts the ambiguities of Ref. [l8| and all 
our results perfectly agree with the magnetization mea- 
surements 2 ^. We also notice that in the real magnetic 
structure deduced from our analysis, there is a small de- 
viation of the Tb magnetic moments from the local Ising 
(111) axes, which is discussed below. 



IV. MONOPOLE PICTURE 

The high- field magnetic structure of Tb2Ti207 can be 
viewed as tetrahedra layers of "1-in, 3-out" type alternat- 
ing with layers of opposite type "3-in, 1-out" . Such local 
configurations are exactly the defects with respect to the 
ice rules, namely the monopoles, which are the building 
bricks of the magnetic monopole representation in spin 
ices^. Thus, the Tb2Ti2 07 magnetic structure consists 
of a 3D array of alternating double layers of monopoles 
and antimonopoles (see Fig. [5]), corresponding to positive 
and negative magnetic charges, and stacked perpendicu- 
lar to [001]. Note that this periodicity is different from 
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FIG. 2. (Color online) Neutron scattering of Tb 2 Ti 2 07 in OT (left panel), Tb 2 Ti 2 7 in 5T (middle panel) and Ho 2 Ti 2 7 in 
2T (right panel) at 1.6 K with H || [110]. Intensity is given in a logarithmic scale. Selected magnetic reflections are highlighted. 
In zero field, only the peaks of the crystal structure are seen. 
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FIG. 3. (Color online) A typical fit for Tb 2 Ti 2 7 at 1.6 K 
and 7T. The calculated integrated intensities (-F<?aic) are plot- 
ted against the experimental ones (F^ hs ). Left panel: Both 
magnetic k — and nuclear intensities. Right panel: Pure 
magnetic intensities of the k — (0, 0, 1) structure. 




FIG. 4. (Color online) Magnetic structures of Tb 2 Ti 2 7 (left 
panel) and Ho 2 Ti 2 7 (right panel) in a [110] field. Blue and 
red balls correspond to opposite magnetic charges. The orien- 
tations of the P moments inside a chain are AFM in Tb 2 Ti 2 7 
and FM-like in Ho 2 Ti 2 7 . 



that of the staggered monopoles structure 2 ^ observed in 
spin ices when the field is along the [111] direction. 

In contrast, in the spin ice Ho2Ti2 7 studied in the 
same conditions, the highest intensity peaks measured in 



neutron diffraction have the lowest intensity in Tb2Ti20 7 
and vice versa (Fig. [2j right panel). The same symmetry 
analysis shows that in Ho2Ti20 7 the j3 chains also order 
antiferromagnetically with respect to each other above 
H « 0.1 T. But here the orientation of the /3 moments 
alternating "in" and "out" combined with the a moments 
alternating in the same way preserve the ice rules in this 
single state (Fig. 2J right panel) 12 i 13 . The refined mag- 
netic structure is close to the schematic X structure 22 , 
with the difference that the moments of the j3 chains are 
lower (8/ie) than the saturation value of 10 /iB 19 - We 
recall that the X structure involves a and /? moments 
of equal magnitudes aligned along the local (111) axes, 
but alternating to keep a net ferromagnetic moment in a 
given chain (FM-like behavior). In this single state, the 
antiferromagnetic coupling of the /? chains together with 
the ferromagnetic coupling of the a chains preserves the 
ice rules. The reduced moment value obtained in the ex- 
periments and the Lorenzian peak shape shown in Fig. [2] 
come from a partial disorder among the f) chains, as dis- 
cussed below. 

Neglecting this effect and using the monopole descrip- 
tion as above, the ordered state of Ho2Ti20 7 spin ice in 
a magnetic field H || [110] contains no monopoles in the 
ground state, since the monopole charges compensate in 
each tetrahedron. That is, the ground state of Ho2Ti20 7 
in a field is "2-in, 2-out" structure, which represents a 
"vacuum" state from the monopole point of view. We 
underline that the main difference between the field- 
induced magnetic structures of Tb2Ti2 7 and Ho2Ti20 7 
comes from their k = (0,0,1) structures, namely, the 
type of ordering within the j3 chains, which is responsi- 
ble for the overall monopolar /vacuum order. These AFM 
k = (0,0,1) structures are described by different mag- 
netic space groups Inna and Inma for Tb2Ti2 7 and 
Ho2Ti20 7 , respectively, whereas both compounds have 
the same magnetic symmetry I mm' a' of their FM-like 
k = structures. 

Starting from these two ordered states mentioned 
above for both Tb2Ti20 7 and Ho2Ti20 7 , elementary ex- 
citations can be viewed as follows. In Tb2Ti20 7 , an exci- 
tation in the double-layered monopolar state is obtained 



FIG. 5. (Color online) Left panel: Double- layered monopolar structure of Tb2Ti2 07 with vacuum pair excitations. Right panel: 
Magnetically vacuum state of Ho2Ti2 07 with monopole pair excitation. 



by reversing the orientation of a single magnetic moment 
between a monopole and an antimonopole (but not be- 
tween two monopoles of the same "type" ) to recover the 
spin ice "2-in, 2-out" structure in the tetrahedra involved. 
This corresponds to a recombination of two monopole 
and antimonopole nearest neighbors, which produces a 
pair of magnetically neutral vacancies. Such vacancies 
can move away from each other (Fig. [5j left panel) with 
little energy cost by further flips of adjacent magnetic 
moments provided they remain within the same j3 chain. 
Hence although monopoles order in three dimensions, 
their excitations only propagate in one dimension along 
the j3 chains. In turn, the a moments with anisotropy 
axes close to the field experience a large Zeeman energy 
and hence cannot flip easily. Alternatively, the propaga- 
tion of a monopole pair in the vacuum state of Ho2Ti2 07 
in a [110] field (Fig. [5j right panel) mirrors that found in 
Tb 2 Ti 2 7 . 

In the high- field ordered states of Tb2Ti207 and 
Ho2Ti207 with H || [110], defects (namely monopoles or 
monopole vacancies) can be created through a small field 
misalignment. They can be evidenced in neutron diffrac- 
tion experiments as deviation from the perfect Bragg re- 
flection, such as a peak broadening or diffuse scatter- 
ing. Increasing the field misalignment and the number 
of defects breaks the monopolar / vacuum order in both 
Tb 2 Ti 2 07 and Ho 2 Ti 2 07 compounds! 3 -^, but the mech- 
anisms are different. Due to its peculiar crystal field 
scheme 24 , Tb 3+ has a finite value of the g± component 
of the Lande tensor, g±/g\\ ~ 1/5 in a field of IT 25 , an 
order of magnitude larger than those in Ho or Dy spin 
ices with Ising behavior. This no n- negligible g± leads to 
the appearance of a magnetic component along the field 
in the j3 chains, so that the magnetic moments deviate 
from their local (111) Ising axes (Fig. 2]). Consequently, 
in Tb2Ti207 the double- layered monopolar order disap- 
pears when increasing the field misalignment, into a spin- 
ice-like configuration involving moments of unequal mag- 
nitude a 12 i 13 . In contrast, in Ho2Ti2 07 the deviation of 



the ordered moments with respect to the anisotropy axes 
is negligible. A huge precision in the field alignment is 
necessary to achieve a fully ordered vacuum state (X 
structure). Such a vacuum state breaks down with the 
misalignment and another ordered state which also obeys 
the ice rules, namely the k = structure, is stabilized. 

As shown previously^ 4 ., dipolar interactions are suffi- 
cient to describe the high- field state of Ho2Ti207. In 
dipolar spin ices (DSI) where the dipolar interaction 
dominates over the exchange, the isolated spin chains 
tend to order ferromagnetically. The simplest model that 
captures the essential features of the behavior of the f3 
chains in a magnetic field is the model of non-interacting 
Ising chains with nn exchange and dipolar interactions 2 ^. 
We adapt this model to Tb2Ti2 07, as shown below. 

The magnetic interaction Hamiltonian is taken as 

H = H ex + H dip = -J J2 { i,j) Ji ■ J j 

i T) c>3 Ji-Jj-3(Ji-Ri j )(J j -R i j) 
^ UIX nn 2^i>j \ Rij \3 

Here Rij = Rj—Ri = | R^ \ Rij , Ri is the position vector 
of atom i with total angular momentum «Jj, R nn is the 
distance between nn atoms, V = /io(#JMB) 2 (47ri? 3 n ) _1 
is the dipolar coupling, gj is the Lande factor, and J is 
the nn exchange coupling with J < for AFM. 

In order to estimate the strength of the dipolar and 
exchange interactions of a single chain in high field (above 
4T), we use our experimental values of the Tb moment 
orientation, which is governed by the crystal field. The 
Tb-/3 moments of the k = (0, 0, 1) structure deviate from 
the local (111) axes by about 13°. In this case, the nn 
exchange interaction is simplified to 

where \i is the Tb magnetic moment. The long-range 
dipolar interactions can be calculated by considering the 
interaction between the atom and its odd neighbors at a 
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FIG. 6. (Color online) Left panel: The structural distortion 
in Tb2Ti2 07 compatible with the observed high- field mag- 
netic structure. The directions of displacements of the axial 
oxygens are shown by white arrows. Right panel: The undis- 
torted structure of Ho2Ti2 07- 

distance i?2m+i = (2ra + 1) R nn or its even neighbors at 
a distance R^m = 2ra R nn as 

H dip * V^(gjfi B )- 2 [l.8 ^(2m + l)" 3 - 1.5 ^(2m)- 3 

« lAA p V^(gjfi B )- 2 « 1.7VfjL 2 (gjfji B )- 2 , 

with Apery constant A p = ^ 1/m 3 w 1.2. 

Similar to dipolar spin ices 2 , to take exchange and 
dipolar interactions into account in Tb2Ti 2 07, we define 
an effective energy scale, H e ff, as 

Heff = (0.7 J + 1.7V) Li 2 (gjii B )- 2 , 

yielding an effective exchange J7 e ff = 0.7 J + 1.7X>. 

So, J7 e ff is expected to be negative (antiferromagnetic) 
when J/V < -1.7/0.7 w -2.4. We estimate the value of 
considering the exchange parameters for Tb2Ti2 07 
reported in the literature. Taking into account the exper- 
imental value of R nn « 3.59 A for Tb2Ti 2 07, the dipo- 
lar coupling is estimated as V « 0.0315 K 9 . There are 
two different published values of J , namely, — 0.083 K 24 
and -0.167K 9 , giving J /V of about -2.6 and -5.3, re- 
spectively. Nevertheless, both of them yield a negative 
J e R and are in agreement with an AFM order within the 
P chain. Note that J « -0.083 K 2 ^ is expected to be 
more reliable as it takes into account the full crystal field 
scheme of Tb 2 Ti 2 7 . 

Thus, the comparison between the strength (and signs) 
of the dipolar and exchange interactions is consistent 
with the AFM order inside the /3 chains. However, the 
model applied in here deals with a single chain and thus 
it cannot explain the full 3D magnetic order. Therefore, 
this oversimplified model should be extended by taking 
into account both the anisotropic exchange and the full 
crystal field scheme. Moreover, numerical Monte Carlo 
studies 14 show that in the X structure of the Ho2Ti2 07 
spin ice, both the FM moment components along a j3 
chain and the AFM orientation of nn chains arise from 
the dipolar interaction. This explanation does not hold 
for Tb2Ti2 07 where an additional mechanism should be 
at play. We propose to consider the influence of a dis- 
tortion, of Jahn- Teller or spin-Peierls type. Structural 



fluctuations were observed in Tb2Ti207 at low tempera- 
ture 27 , possibly precursor of a T = Jahn- Teller transi- 
tion. A spontaneous Jahn- Teller distortion could be the 
main ingredient to explain the spin-liquid ground state of 
Tb2Ti207 in zero fieldii. Exchange driven spin-Peierls 
distortions have also been considered as possible mecha- 
nisms to relieve the degeneracy of the pyrochlore lattice, 
in spinels compounds especiall y 28 ! 29 . 

Up to now, the structural distortion in Tb2Ti207, ei- 
ther spontaneous 27 , field induced 30 or induced by pres- 
sure/stress^ 1 -, has not been characterized. In our case 
when the field applied along the [110] axis its symmetry 
should reflect that of the magnetic order. We derived the 
distortions compatible with the observed magnetic modes 
by symmetry analysis. Interestingly, we found that only 
one distorted structure has a symmetry compatible with 
that of the magnetic order of Tb2Ti207 in high field, and 
that this structure involves displacements of the oxygen 
ions exclusively. We recall that in the pyrochlore struc- 
ture, each rare-earth ion is surrounded by eight oxygens, 
two "axial" ones along the local (111) axis and six "pla- 
nar" ones lying in a buckled perpendicular plane. The 
symmetry constraints on the allowed displacive modes 
are very different for axial and planar oxygens. The dis- 
placement of axial oxygen is allowed only along the [001] 
axis and occurs in such a way that half of the oxygens 
move parallel and half antiparallel to this axis (see Fig. [6j 
left panel). The displacive mode of planar oxygens is 
more complex but can be also explicited. However, we 
believe that the displacements of planar oxygens have 
less effect on the magnetic interactions and on the crys- 
tal field of rare-earth pyrochlores, since these oxygens are 
typically about 13 % farther from rare-earth ions than the 
axial ones. Hence, we suggest that in the monopole pic- 
ture proposed here for Tb2Ti2 07 in high magnetic field 
along [110], oxygen displacement should result from the 
attraction (repulsion) of monopoles with opposite (iden- 
tical) magnetic charges. For comparison, we also show 
the reference structure of Ho2Ti2 07 (see Fig. [6l right 
panel) where such distortions are expected to be negligi- 
bly small, due to its different crystal field scheme. 



V. CONCLUSION 

In conclusion, we show that under a high magnetic 
field applied along a [110] direction, the magnetic struc- 
ture of Tb2Ti2C>7 can be viewed as a 3D arrangement 
of monopole and antimonopole double layers, in contrast 
to Ho2Ti207 where the high- field magnetic structure be- 
haves as a vacuum (monopole- free) state. A toy model for 
an isolated chain shows that in Tb2Ti207 the exchange 
interaction overcomes the dipolar one, which is consistent 
with AFM order in the /3 chains, but understanding the 
full 3D magnetic structure of Tb2Ti207 requires to take 
into account the real crystal field, anisotropic exchange, 
dipolar interactions and the influence of a distortion. Fi- 
nally, symmetry analysis allows us to derive a distortion 
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model compatible with the observed magnetic structure of Tb2Ti 2 07 which can be related to the appearance of 

the double-layered monopolar state. 
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